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The strange and light quark contributions to the nucleon mass from Lattice QCD 
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We determine the strangeness and light quark fractions of the nucleon mass by computing the 
quark line connected and disconnected contributions to the matrix elements mq{N\qq\N) in lattice 
QCD, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We 
simulate np = 2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing 
a ~ 0.073 fm. The renormalization of the matrix elements involves mixing between contributions 
from different quark flavours. The pion-nucleon a"-term is extrapolated to physical quark masses 
exploiting the sea quark mass dependence of the nucleon mass. We obtain the renormalized values 
CttN = (38 ± 12) MeV at the physical point and fr, = cTs/m^ — 0.012(14)^3° for the strangeness 
contribution at our larger than physical sea quark mass. 

PACS numbers: 12.38.Gc,13.85.-t,14.20.Dh 
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I. INTRODUCTION 

Most of the nucleon's mass is generated by the sponta- 
neous breaking of chiral symmetry and only a small part 
can be attributed directly to the masses of its valence 
and sea quarks. The quantities 



/t, ^mq{N\qq\N)lmn 



(1) 



parameterize the fractions of the nucleon mass ttin that 
are carried by quarks of flavour q. Almost all visible mat- 
ter of the universe is composed of nuclcons and yet little 
is known experimentally about these quark contributions 
to their mass. 

The scalar matrix elements 



mq{N\qq\N) = m^fj 



(2) 



also determine the coupling strength of the Standard 
Model (SM) Higgs boson (or of any similar scalar parti- 
cle) at zero recoil to the nucleon. This then might couple 
to heavy particles that could be discovered in LHC ex- 
periments, some of which are dark matter candidates [l| . 
The combination m^J^qfT,, 9 G d, s}, will appear 
quadratically in this cross section that is proportional to 
I/nP where 



777, js^ ^— ^ m . /.I ^— ^ m . 
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with the couplings aq oc mqjvny/. Due to the trace 
anomaly of the energy momentum tensor one obtains 



ha 



(4) 



Note that the coupling /n of Eq. ([3]) only mildly depends 
on the masses of heavy quark flavours [2, Q • 

The fjg-terms are also needed for precision measure- 
ments of SM parameters in collisions at the LHC. For 
instance the resolution of a (hypothetical) mass differ- 
ence between the and W~ bosons is limited by our 
knowledge of the asymmetries between the up and down 
as well as between the strange and charm sea quark con- 
tents of the proton Q . An accurate calculation of these 
quantities will help to increase the precision of SM phe- 
nomenology and to shed light on non-SM processes. 

The light quark contribution, the pion-nucleon cr-term, 
is defined as 



= 



9mN 



drud 



(5) 



where mps denotes the pseudoscalar mass. Some time 
ago, employing dispersive analyses of pion-nucleon scat- 
tering data, the values [sj g-^n = 45(8) MeV and Q 
CttN = 64(7) MeV were obtained while a calculation in 
the framework of Oip^) heavy baryon chiral perturbation 
theory resulted in Q ct^n = 48(10) MeV. A recent covari- 
ant baryon chiral perturbation theory (B^FT) analysis 



of the available pion-nucleon scattering data fs'l resulted 
in the value ctttN = 59(7) MeV. Even less is known about 
the strangeness contribution (Ts- Since no elastic Higgs- 
nucleon scattering experiments exist all phenomenolog- 
ical estimates largely depend on modelling. Therefore 
input from lattice simulations is urgently required. 

We have witnessed an upsurge of interest in calculat- 
ing flavour singlet matrix elements recently, either di- 
rectly , by calculating the corresponding quark line 
connected and disconnected terras, or indirectly, via the 
Feynman-Hellmann theorem |13l - [l9j . High statistics sim- 
ulations including light sea quarks mean that reasonable 
signals can be obtained for disconnected terms. Simi- 
larly, the small statistical uncertainty on baryon masses 
as functions of the quark masses enable reasonable fits to 
be made. Ideally, the results of both approaches should 
agree. 

Preliminary results of this study were presented at past 
Lattice conferences [l^, [2l| . This article is organised as 
follows. In Sec. |lT] we detail the gauge configurations, 
simulation parameters and methods used. In Sec. IIIII 
we then explain how the lattice results are renormalized 
and finally we present our results in Sec. llVi before we 
summarize. 



II. SIMULATION DETAILS AND METHODS 

We simulate np = 2 non-perturbatively improved 
Sheikholeslami-Wohlert fermions, using the Wilson gauge 
action, a.i (3 — 5.29 and n — Hud — 0.13632. Setting the 
scale from the chirally extrapolated nuclcon mass, we ob- 
tain {221 the value tq = 0.508(13) fm for the Sommer 
scale, in the physical limit. This results in the lattice 
spacing 

= (6.983 ± 0.049) r^^^ = (2.71 ± 0.02 ± 0.07) GeV , 

(6) 

where the errors are statistical and from the scale setting, 
respectively. The r^/a ratio is obtained by chirally ex- 
trapolating the QCDSF /3 = 5.29 simulation points 23). 
An extrapolation of the axial Takahashi-Ward identity 
(AWI) mass yields the critical hopping parameter value 

Kcsoa = 0.1364396(84) . (7) 

In addition to K^d = i^sca — 0.13632, we realize the 
valence K-values, — 0.13609 and Kg = 0.13550. The 
corresponding three pseudoscalar masses read 

mps.ud = (0.1050 ± 0.0003) 

= (285 ± 3 ± 7) MeV , (8) 

mps,m = (449±3±ll)MeV, (9) 

mps,s = (720±5±18)MeV. (10) 

The strange quark mass was fixed so that the above 
value for mps,s is close to the mass of a hypothetical 
strange-antistrange pseudoscalar meson: (m^± -I- 
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FIG. 1. Quark line connected (top) and disconnected (bot- 
tom) three-point functions. We have omitted the relative mi- 
nus sign between the diagrams. Note that for scalar matrix 
elements, the vacuum expectation value of the current inser- 
tion needs to be subtracted [qq ^ qq — {qq}), see Eq. (|12|l . 

m2±)i/2 ~ 686.9 MeV. We investigate volumes of 32^64 
and 40'^64 lattice points, i.e., Lmps — 3.36 and 4.20, 
respectively, where the largest spatial lattice extent is 
L ss 2.91 fm. We analyse 2024 thermalized trajectories 
on each of the volumes. To effectively eliminate auto- 
correlations for the observables that we are interested in, 
bin sizes of eight are found to be sufficient. 

The matrix element {N\qq\N) is extracted from the 
ratio of three-point functions, see Fig. [1] to two-point 
functions at zero momentum. Defining 

TrtA = ^TrA(x,t;x,0, (11) 

X 

we can write the disconnected part as 

«-,M) = (T,,(M--.)>-<^^i^a, 

(12) 

where M is the lattice Dirac operator for the quark 
flavour q of the current and C2pt(if) denotes the two- 
point function of the zero momentum projected proton 
connecting the source time ti = with t{. Note that, 
unlike its expectation value, C2pt computed on one con- 
figuration will in general also have an imaginary part. 
This means that we can reduce the variance of the above 
expression by explicitly setting this to zero, using the 
relation ImTr((Af ~^1) = that follows from the 75- 
Hermiticity, = 'j^M'j^. 

In the limit of large times, if ^ i ^ 0, 

R'^'^tf, t) + i?^°"(if , t) {N\qq\N) , (13) 

where the term i?™"(tf, t) = C3pt(tf,t)/C2pt(tf) only con- 
tributes for q € {u,d}. Quark field smearing (see below) 
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FIG. 2. Dependence of W^^" on ti for smeared-smeared (SS) 
and smeared-point (SP) two-point functions, together witii 
the fit result. 

at the source and the sink significantly enhances the cou- 
pling of the nucleon creation and destruction operators 
to the physical ground state. Still, the time distances 
between the source and the current insertion t as well as 
between the current and the sink t{ — t need to be taken 
sufficiently large to suppress excited state contributions. 

We find the nucleon smeared-smeared effective masses 
to be constant for t > 8a. It suffices if excited state 
effects are much smaller than the statistical errors. These 
errors however are expected to be substantially larger 
for the disconnected three-point function than for the 
nucleon two-point function. Thus, we set the time of the 
current insertion to a smaller value t = 4a ~ 0.29 fm. 
The method that we apply requires us to fix t, but t{ can 
be varied. If the i?'^"^ data were not constant for tf >2t — 
8a then this would have implied that our choice of t was 
too ambitious. Fortunately, employing sink and source 
smearing, we find the asymptotic limit to be effectively 
reached for t{ > 5a and compute the matrix elements by 
fitting the above ratios for > 6a « 0.44 fm to constants. 

As an example, in Fig. [5] we display the disconnected 
ratio for strange valence and current quark masses as a 
function of tf for smeared-smeared three- over two-point 
functions for 40'^64 lattices, together with this fit result. 
In addition we show the corresponding smeared-point ra- 
tio that converges towards the same value, giving us ad- 
ditional confidence that t was chosen sufficiently large 
to warrant ground state dominance within our statistical 
errors. 

Based on Ref. [13] we improve the overlap of our nu- 
cleon creation operator with the ground state by applying 
Wuppertal-smearing [2^ 

= (4-^' + ^ E u.Ar:^ (14) 

to quark fields (/>, where we set S = 0.25 and use 400 
iterations. We replace the spatial links Ux.j above by 



APE-smeared ^ finks 

7-r(n) _ p l^rr(""l) J_ 7-r(«-l)rr("-l) 7-r(n-l)t 1 

V / 

(15) 

where i S {l,2,3},j € {±1,±2,±3}. Psu(3) denotes a 
projection operator into the SU(3) group and the sum is 
over the four spatial "staples", surrounding Ux,i- We em- 
ploy 25 such gauge covariant smearing iterations and use 
the weight factor a = 2.5. For the projector we somewhat 
deviate from Ref. [l^ and maximize ReTr [A^Pgu(3)(A)], 
iterating over SU(2) subgroups. The connected part, for 
which the statistical accuracy is less of an issue, is ob- 
tained with a less effective smearing at the larger, fixed 
value t[ = 15a, varying t. 

We stochastically estimate TrtM~^. For this pur- 
pose we employ N complex Z2 noise vectors, \r]^)t, 
i = 1^...^N, whose spacetime (Ei spin (8) colour com- 
ponents e'"^ carry uncorrelated random phases (/) G 
{±7r/4, ±37r/4} at the time t and are set to zero else- 
where, to reduce the noise (partitioning [13] )■ 

Solving the linear systems 

M\s')t = W)t (16) 
for |s')f we can then substitute, 

1 ^ / 1 \ 

Tr^^^s' = ^ E = Tr,M-^ + O i^-^j . 

(17) 

The inner product is only taken over three-space, spin 
and colour indices. In the case of the scalar matrix el- 
ement it is relatively easy to push the stochastic error 
below the level of the inherent error from fluctuations 
between gauge configuration^ [2^. Therefore, here we 
do not need to employ the Truncated Solver Method 
(TSM) [2^ and do not exploit the hopping parameter 
expansion. Instead, to reduce the dominant gauge er- 
ror, we compute the nucleon two-point functions for four 
equidistant source times on each gauge configuration. 
In addition we exploit backwardly propagating nucle- 
ons, replacing the positive parity projector ^(1 + 74) by 
^(1 — 74) within the nucleon two-point function, C2pt- 
Consequently, the noise vectors are seeded on 8 time 
slices simultaneously, reducing the degree of time par- 
titioning. We find this not to have any adverse effect on 
the stochastic error. In addition to the 48 (4 timeslices 
times 4 spinor components times 3 colours) point-to-all 
sources necessary to compute the two-point functions we 
solve for TV = 50 noise vectors per configuration and cur- 
rent quark mass. 



Note that both error sources will scale in proportion to 1/ y/n^onl- 
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III. RENORMALIZATION 

In the continuum, for light quark flavours q, the CTg- 
terms are invariant under renormalization group trans- 
formations. However, Wilson fermions explicitly break 
chiral symmetry and this enables mixing not only with 
gluonic contributions but also with other quark flavours. 
Note that due to the use of a quenched strange quark, 
the renormalization of the corresponding matrix element 
is particularly large and needs to be studied carefully. A 
consistent 0{a) improvement of the quark scalar matrix 
elements requires the inclusion of the gluonic operator 
aGG. We have not measured this as yet. Therefore we 
will neither include any 0{a) improvement of the renor- 
malization constants nor of the scalar current. However, 
we will account for the mixing between quark flavours. 

We follow the procedure outlined in Ref. [2§|, see also 
Sec. 6 of Ref. |30|. The same result can be obtained 
by taking derivatives of the nucleon mass with respect 
to the sea quark masses [3ll | via the Feynman-Hellmann 
theorem. This also holds for the case of a partially 
quenched strange quark [isj . In the renormalization 
the strangeness matrix element will receive subtractions 
not only from light quark disconnected but also from the 
numerically larger connected diagrams. This has first 
been pointed out in Ref. [s^ . 

For Wilson actions the vector Takahashi-Ward identity 
(VWI) lattice quark mass is given by 



1 / 1 

2a 



1 



c,sca 



In the partially quenched theory one can define a 
Kc,vai(/3, Ksca) 7^ 'tc,soa(/3) as the K-valuc at which the va- 
lence pseudoscalar mass vanishes. 

We distinguish between singlet and non-singlet quark 
masses that will renormalizc differently. In the case of the 
theory with np — 2 + 1 {rriu — md — mud) sea quarks, 
we have a mass term 



CO 0.04 - 



c 




0.01 0.02_ 0.03 0.04 
ma 

FIG. 3. Determination of the slope Z^/Z"!^ — 1 + az, see 
Eqs. ([25 )) and (p7)l . 



In the partially quenched theory with mass degenerate 
sea quarks Eq. 1^1^ results in 

m^^^^i) = z:Mm, (22) 

m-"(A.) - = (m)(^ - ™™') , (23) 



(18) where we introduce a VWI valence quark mass through 



TO 



val 



1 



f. 



1 



1 



At Kval — Kc.val: "m, 

Zns 



2cL \ /^val ^^c,i 

vai,ren yauishes SO that we obtain 

val 



(24) 



TO 



TO — TO 



TO 



TO 



Lrn = mud{uu + dd) + nisSS 

TO = i(2TO„d + TOs) , 



V.val 



(19) 
(20) 



where ip — {u,d,s); the lattice singlet quark mass to is 
given by the average of the sea quark masses. 

A renormalized quark mass for the flavour g at a scale 
H is given by 

TO-"(m) = ZM^+Z:^'i^^)im,-^) : (21) 



/^c.sca 



= : 1 + az . 



(25) 



The non-singlet mass above can also be obtained from 
the bare AWI mass. 



= Z^to"^ = 



(26) 



that renormalizes with the ratio of the renormalization 
constants over Zp" of the non-singlet axial and pseu- 
doscalar currents. This results in the alternative prescrip- 
tion to Eq. ((25)) . 



where singlet quark masses renormalize with a renor- 
malization constant Z^(/i) and non-singlet combinations 
with Z"j'*(^). Notice that the ratio Zf„/Z"j'* is indepen- 
dent of the scale /i but depends on the lattice spacing a, 
through the coupling Q;s(a~^). 



7ns 



2a TO 



AWI 



yns yris — i 
^ P /^sca 



(27) 



where is calculated at Ksca- The required com- 

bination of scalar, pseudoscalar and axial non-isosinglet 
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renormalization factors for our simulation with np = 2 
sea quarks at (3 = 5.29 reads (s^l 

^r/(^m ^p ) = 0-988 ± 0.031 . (28) 

The two methods that differ by terms of 0{a) are illus- 
trated in Fig. [3] for our configurations with np ~ 2 mass 
degenerate sea quark flavours. 
The results read 

_ r 0.411(13) , VWI 

^ \ 0.369(22) , AWI. ^ 



Corrections to these values at a non-zero quark mass m 
will be of 0{ma); the differences between the two defini- 
tions are indicative of 0{a) effects. We remark that az 
is of leading order in the strong coupling parameter 
and it will therefore decrease with the lattice spacing a. 
For ?iF = 2 + f sea quarks Eq. (PT|) amounts to 

m-"(Ai) = Zr(M)"^, + (KM - ^r(M))m (30) 
= (a*) (™<? + Oizm) , 

where az is defined in Eq. (P5|) . This can be written as 



m„(^) 



-z;-(M,a) 



V 



az(a) 
^ 3 
az(a) 
3 

az(a) 
3 



az(a) 
3 

^ az(a) 




m„(a) 
md{a) 
ms{a) 



lat 



(31) 



where the lattice quark masses on the right hand side are defined by the VWI, Eq. p^. For clarity we have included 
the lattice spacing dependence above, which below we will drop again. The sum rrf™{^){N\qq\Ny°'^{ii) is invariant 
under renormalization group transformation^ Therefore the scalar lattice matrix elements will renormalize with the 
inverse matrix above: the different quark contributions will mix in the mass non-degenerate case. 

We now turn to the situation of interest of rtp = 2 light sea quarks, with a quenched strange quark. This means 
that the singlet mass m'*** — (mjf' -I- mJj'*)/2 will not depend on the strange quark mass anymore. (The superscript 
"lat" has been added for clarity.) Eq. (150)) can now be written as 



ms{ii) 



QZ. 

2 

az 
2 



ctz 
2 



ctz 
2 



az 
2 



oz. 
2 





/ \ lat 




•) 


■md 




[ ms / 



(32) 



where the presence of the strange quark is not felt by the sea quarks. However, the definition of involves Kc 
Inverting the above matrix yield^ 



/ {uu) {p) \ . 
(rfrf)(M) =77—^ 



1 



I Qz. Oz. 

2 2 
Qz 1 I Qz 
2 -L -p o 
ctz 
2 



For the light quark matrix element (i.e. the (iTrN-term) this means that 
mr(M) + "^r(M) 






^ l + az 




(33) 



J at 



-{N\uu + dd\NY'"^(^i) = 



^lat 



{N\uu + dd\N) 



lat 



(34) 



while for the strangeness matrix element we obtain 

K(7V|ss|7V)]"" = 



lat 



lat 



lat\ 
d ) 



{N\ss\N) 



lat 



az 



2(1 + az) 



(iV|uu-|-d(i|7V) 



lat 



(35) 



Note that these scalar matrix elements are differences of scalar 
currents qq within the nucleon, relative to their vacuum expecta- 
tion values. Therefore, unlike the chiral condensates alone, they 
do not undergo any additive renormalization and do not mix 
with an a~^l term. 

Here and occasionally below we omit specifying the external state 
(in our case \N)) in cases where identities between hadronic ma- 
trix elements are independent of this state. 



Again, the lattice strange quark mass is defined as m 



lat 



5^)/ (2a). The same renormalization pattern 



can also be derived, employing the Feynman-Hellmann 
theorem 



It is evident from Eq. (|33p that the so-called y-ratio 
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mps/GeV^ 



FIG. 4. The unrenormalized mass fraction f^^, as a function 
of the current pseudoscalar mass on the L = 32a lattices, at 
the smallest valence mass mps ~ 285 MeV. 



renormalizes as follows, 

2{N\ss\NY 



y 



{N\uu + ddlNY'"^ 

2(A^|ss|7V)i^' 



(1 + az) 



{N\uu + dd\N) 



lat 



az 



(36) 



IV. RESULTS 

As discussed in Sec. [TT] above we employ three hopping 
parameter values, Kud — ^>■scs^ — 

0.13632, ^ 0.13609 
and Kg = 0.13550, that correspond to the pseudoscalar 
masses tops w 285 MeV, 450 MeV and 720 MeV, respec- 
tively, see Eqs. (O - (|TU)) . We use all these k- values for 
the valence quarks (^vai) as well as for the current inser- 
tions qq (ftcur)- This amounts to nine combinations for 
the disconnected ratios R'^^^ while for the connected part 

We will explore the dependence of the lattice matrix 
elements on the current quark mass and investigate finite 
size effects, using a partial summation method. Subse- 
quently, unrenormalized and renormalized valence and 
sea quark contributions are studied. Finally, we com- 
pute the light cr-term, the mass contributions fx, and 
/tg and the y-ratio. 




FIG. 5. Eq. H12|l . partially summed up to a maximum spatial 
distance from the source Xmax, see Eq. (P7)l . 



contribution in the hopping parameter expansion is pro- 
portional to the plaquette. To confirm this expected sat- 
uration we compute the scalar matrix element for three 
additional current quark masses up to and above the 
charm quark mass TOc on a subset of 576 32^64 lattices, 
using our smallest valence quark mass. We remark that 
for an 0{a) improvement of the operator qq, mixing with 
aGG needs to be considered. These improvement terms 
will become large if to^o is big. In fact the renormalized 
becomes negative at about twice the strange quark 
mass if we neglect such effects. We have not investigated 
gluonic contributions as yet and therefore the behaviour 
at large masses is beyond the scope of the present study. 
For very small quark masses we would expect /t, oc to^. 
Indeed, these expectations are confirmed for the unrenor- 
malized /j?', as can be seen in Fig. 21 Note that the third 
data point from the left corresponds to the strange quark 
mass while the charm quark can be found at the value 



'"•PS 
eye. 



.9GeV^ Thearctan -curve is drawn to guide the 



B. Finite size effects: partial sums 

In order to investigate finite size effects we find it 
worthwhile to replace the numerator of Eq. (IT^ by a 
partial sum: 



E 

|X|<2;„„: 



(C2pt(tf)Tr[M-i(x,t;x,0])^ 

(C^2pt(tf)) 



(37) 



A. Dependence of the lattice matrix elements on 
the current quark mass 

From the heavy quark expansion it is evident that 
Jt^ oc {N\GG\N)/mj<! oc Jtg for ruq — >■ oo. This can also 
be seen on the lattice where the leading non-vanishing 



The subscript c (connected part) indicates that we sub- 
tract the product of the two individual vacuum expecta- 
tion values from the numerator. Since no zero momentum 
projection is performed at the source of the two-point 
function, that resides at the spatial position x = 0, the 
result will depend on the cut-off Smax- We expect the 
summand at large |x| to fall off exponentially in |x| so 
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TABLE I. The disconnected contribution to the scalar lattice 
matrix elements for diflterent K-combinations. 







V 


{N\qq\Nfi^, 




0.13550 


32^64 


2.01(21) 






40^64 


2.17(25) 


0.13550 


0.13609 


32^64 


2.27(22) 






40^64 


2.43(27) 




0.13632 


32^64 


2.38(23) 






40^64 


2.55(29) 




0.13550 


32^64 


1.97(20) 






40^64 


2.06(25) 


0.13609 


0.13609 


32^64 


2.19(22) 






40^64 


2.28(26) 




0.13632 


32^64 


2.21(23) 






40^64 


2.36(28) 




0.13550 


32^64 


1.96(23) 






40^64 


1.93(27) 


0.13632 


0.13609 


32^64 


2.06(24) 






40^64 


2.05(29) 




0.13632 


32^64 


1.67(26) 






40^64 


1.86(31) 



that these values will eventually not contribute to the 
signal anymore but just increase the statistical noise. At 
very large spatial volumes one may therefore consider to 
perform such a partial sum only, thereby reducing the 
statistical error, and to estimate the induced bias by pa- 
rameterizing the asymptotic fall-off. 

We display the partial sums for the L — 32a and L = 
40a lattices for the strangeness current at Kvai = '*sca 
in Fig. [5] We do not detect any statistically significant 
dependence of the curves on the value of tf and show 
the results obtained at U — 6a. At small a^max we see 
the naively expected x^^^^^^ volume scaling. This becomes 
flatter around Xmax ~ 8a but only saturates to a constant 
once the boundary of the L = 32a box is hit at Xmax = 
16a. Increasing Xmax beyond this value means that in 
the case of this smaller lattice only the lattice "corners" 
are summed up. However, the L = 40a data saturate 
at the same distance, rather than at 20a, indicating that 
indeed the nucleon is well accommodated within this box 
size and that finite size effects are small. 



C. Sea and valence quark contributions 

The results of the bare disconnected scalar matrix el- 
ements for the two volumes (and three quark masses 
niq < TTLs) are displayed in Table ID For disconnected 
terms Kcur can differ not only from Kgoa but also from the 
Kvai of the nucleon's valence quarks. In Fig. |6]we display 
the dependence of the unrenormalized /j?' -values on the 




0.08 - 

0.07 ' 1 1 ' ' 1 1 

mps=m^ 1 2 3 

(mpsro) 

FIG. 6. The unrenormalized strange quark mass fraction /;^*, 
as a function of the valence pseudoscalar mass, for the two 
volumes. 

valence quark mass of the proton for both volumes, to- 
gether with linear chiral extrapolations. The right-most 
data points correspond to the strange quark mass and 
the left-most points to the present sea quark mass. The 
volume dependence is not significant and neither are the 
differences between the values obtained at the smallest 
mass point and the chirally extrapolated numbers. The 
results need to be renormalized and this is possible at 

'^val — /^soa — 0.13632. 

To enable the calculation of the light cr-term and 
the renormalization of the strangeness matrix element, 
we also compute the connected contribution for Kcur = 
Kvai = f^sca = 0.13632, using the traditional_ sequential 
propagator method. We obtain {N\uu -f- dd\NY^^^ ~ 
8.43(73) and 8.35(43) for the L = 32a and L = 40a lat- 
tices, respectively. This means that at the pseudoscalar 
mass mps ~ 285 MeV the relative contribution of the 
disconnected matrix element reads 

^lat ^ {N\uu + dd\Ny^^ ^ { 0.284(36) , L = 32a 
{N\uu + (idliV)!^' I 0.308(37) , L = 40a . 

(38) 

The unrenormalized values of Table U seem to be fairly 
independent of the current quark mass. The ratio of 
the strangeness matrix element over a light sea quark 
contribution undergoes the renormalization 

2{ssy-- ^ 2{-ssy-'~j^{uu + ddy-' 

{uu -f ddy^^ {uu + Jd)[fi* - {uu + ddy^^ ' ^ ' 

The renormalization of the numerator is obtained from 
Eq. p3|) while the denominator can be split into two 
parts that renormalize with and with , respec- 
tively: {uu + dd) — {uu + dd) con- The SU(3)f flavour 
symmetry of the unrenormalized sea obviously cannot 
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disappear when subtracting the same (large) terms from 
the numerator and the denominator. However, this sub- 
traction resuhs in large statistical uncertainties. For in- 
stance on the large volume the above ratio reads 1.7 ±5.5 
with an additional systematic uncertainty of 0.5 from the 
value of the renormalization parameter az, Eq. ([2^ . We 
conclude that the renormalized sea is SU(3)f symmetric 
within a factor of about five. 

The disconnected fraction r of Eq. p8D will undergo 
the renormalization 



m^[MeV] 



r'™ = (1 + az)r'^' - az = 0.024(5)+ 



+29 
9 • 



(40) 



The value quoted is obtained on the L = 40a volume us- 
ing the VWI prescription, with a systematic error that 
incorporates the difference between the two determina- 
tions of the renormalization constant ratios Eq. and 
their respective uncertainties. For the renormalized y- 
parameter that is defined in Eq. p6p this implies that 



2(ss)"' 



y 



{uu + dd) 



rcn 
dis 



(41) 



where the approximation holds within a factor of five, see 
also Eq. (gT]) below. 



D. The light and strange cr-terms 

Combining all information results in the renormalized 
values for the pion-nucleon a-term at the simulated sea 
quark mass 



CPSN 



_ ( 0.0378(39)a-i = 0.264(26)(2)r 
~ I 0.0389(36)a-i = 0.272(25)(2)t 



L = 32a 
L = 40a 
(42) 



1.5 



1.0 



0.5 - 



0.0 - 




0.2 



0.4 

(romps)^ 



0.6 



FIG. 7. Extrapolation of crpsN/fips to the physical point [221 
using covariant BxPT. The open symbols correspond to the 
values that we directly obtain at mpsro ~ 0.73 on the L = 40a 
volume (left) and for L = 32a (right). The broad error band 
is obtained when ignoring this constraint. The horizontal line 
denotes the (constant) leading order expectation. 



TABLE II. QCDSF pseudoscalar and nucleon masses [23l at 
/3 = 5.29 (a-i ^ 2.7 GeV) and /3 = 5.40 (a"^ « 3.2 GeV). 



V 



amps 



arriN romps J^oA^n 



5.29 0.13620 24^48 0.1552(6) 0.467(5) 1.084(9) 3.26(4) 

5.29 0.13632 24^48 0.1112(9) 0.425(6) 0.776(8) 2.97(5) 

5.29 0.13632 32^64 0.1070(4) 0.390(5) 0.747(6) 2.72(4) 

5.29 0.13632 40^64 0.1050(3) 0.381(3) 0.733(6) 2.66(3) 

5.40 0.13660 32^64 0.0845(6) 0.353(7) 0.700(8) 2.92(7) 

5.40 0.13660 48^64 0.0797(3) 0.314(5) 0.660(7) 2.60(5) 



at the two different volumes, where the second error is 
due to the uncertainty of the chirally extrapolated tq- 
value. The above cr-term can also be obtained from the 
derivative of the nucleon mass with respect to the log- 
arithm of the light quark mass. Using the fact that at 



small mps, 



CTPSN 



rrid <x TOpg , we can write 



rrid 



drud 



'PS 



drripc^ 



(43) 



To leading order in chiral perturbation theory 
drrif^ / drripg = const. This linear assumption suggests to 
multiply the result Eq. (|42|) with the ratio phys/™PS 

to obtain the physical cr-term crjjjjf^'° = 0.064(6)rQ"^ — 
25(3)(1) MeV. In general, however, higher order correc- 
tions will lead to some curvature. 

Fortunately, we do not only know the cr-term aX k, — 
0.13632 but also the nucleon mass [22!| at other values of 
K^^i = Ksea, at ^ = 5.29 {a-^ « 2.71 GeV) and at ;9 = 5.4 
(a-i « 3.22 GeV), see Table IE A combined 0(p'') co- 
variant B^PT [34!| fit to these data within the window 



250 MeV < mps < 430 MeV, imposing our directly ob- 
tained value of cTpg n as an additional constraint, results 
in the preliminary number [2^ 



^phys 



(38 ± 12) MeV 



(44) 



at the physical point. The error includes both the statis- 
tical uncertainty of the fit and the systematics from vary- 
ing the low energy parameters C2, C3 and ^3 within their 
phenomenologically allowed ranges [S^-H^- ^ detailed 
analysis will be presented in Ref . . We display the re- 
sult of this extrapolation in Fig.[7]for the ratio ctpsn/wps 
in units of tq, together with our direct determinations. 
The broad error band indicates the result of the same fit, 
without using our constraint at mps ~ 285 MeV. 

We now use Eq. (PSI) with az given in Eq. ([^ to ob- 
tain the renormalized strangeness matrix element from 
the values given above. This amounts to subtracting 
numbers of similar sizes from each other. There is no 
noticeable finite size effect between the 32^ and 40'^ vol- 
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umes. For our simulation point at a low pion mass mps ~ 
285 MeV we obtain the values, a[ms(A^|ss|A^) = 
0.005(6) and 0.008(6) for the two determinations of the 
renormalization parameter az from the VWI and AWI, 
respectively. 

Of particular phenomenological interest is the dimen- 
sionless strange quark contribution to the nucleon mass 



Finally, we also predict the gluonic and heavy sea 
quark contribution /y^ of Eq. Q , 



fn = 



[ms{N\ss\N)] 
mN 



= O.G12(14)tJ°. 



(45) 



Again, we quote the value obtained from the VWI pre- 
scription, with a systematic error that incorporates the 
difference between the two determinations of the renor- 
malization constant ratios and their respective uncertain- 
ties. This may be indicative of 0{a) effects. The prob- 
lem of large cancellations cannot be overcome easily. One 
needs to get closer to the continuum limit so that az ap- 
proaches zero. For instance, at /3 = 5.40, az ~ 0.2 [29l |. 
significantly reducing the subtraction of the connected 
diagram (and probably the value of (iV|ss|iV)'''* that will 
contain a smaller light quark contribution). 

The result obtained is interesting insofar as it suggests 
a scalar strangeness of less than 4 % of the nucleon mass, 
as = 121^5 MeV. In spite of the relative enhancement 
by the ratio ms/rriud > 25 this is not bigger than the 
pion-nucleon cr-term above. This is quite consistent with 
the finding of Eq. (HUl) of a tiny renormalized light sea 
quark participation in ctpsn- We remark that taking 
the combination m|f'*(A^|ss|A^)'^' without the proper sub- 
traction would have resulted in Jt, ~ 0.12, even bigger 
than the light quark mass contribution of about 0.09, at 
our light quark mass value that exceeds the physical one 
by a factor of about four. Neglecting the mixing with 
light quarks in the renormalization is probably the main 
reason why this contribution was overestimated in the pi- 
oneering lattice studies, see e.g. Ref. [^7] and references 
therein. Early results are also summarized in Ref. (T9l | 

We can constrain the scale-independent y-ratio of 
Eq. dSlD, 



(1 az) 0.333(36) - az , L 
(1-haz) 0.320(33) -az, L 

0.059(37)(28) , L = 32a 
0.041(37)(29) , L = 40a, 



32a 
40a 



(46) 
(47) 



where the errors are statistical and the difference between 
the two determinations of az , respectively. Again, as the 
central value, we have taken the result from the VWI 
renormalization factor. From our determination of the 
pion-nucleon cr-term we know that the denominator of 
Eq. will increase by a factor 1.4-1.5 when extrap- 
olated to the chiral limit. Based on the weak observed 
dependence of (A^|ss|A^)'''* on the valence quark mass, 
see Fig. [SJ we would expect the numerator to exhibit a 
less pronounced quark mass dependence. Thus a 95 % 
confidence-level upper limit on the y-parameter y < 0.14 
should also apply at physically light sea quark masses. 



/tg = 1 = 0.951 



+20 
-27 ■ 



(48) 



This means that the light and strange quark flavours con- 
tribute a fraction between 3 % and 8 % to the nucleon 



V. SUMMARY 

We directly calculate the light quark and strangeness 
CT-terms on lattices with spatial extents up to Lmps « 
4.2, a lattice spacing « 2.71 GeV and a pseu- 
doscalar mass mps « 285 MeV. At this mass point and 
lattice spacing the quark line disconnected contribution 
amounts to a fraction of r'°* w 0.3 of the full result. Af- 
ter renormalization however we find this number to drop 
below the 5 % level, see Eq. pOt . 

At our fixed mass point we obtain the renormalized 
values, (Ts = 12l^^ g^^^ ^^^^ ^ 106(11)(3) MeV, 
for the strangeness and light quark cr-terms. Assuming 
the latter value to depend linearly on mpg, as predicted 
by leading order chiral perturbation theory, this corre- 
sponds to 25(3) MeV at the physical point. However, 
from nucleon mass data obtained at pseudoscalar masses, 
250 MeV < mps < 430 MeV, it is clear that there exists 
a non-vanishing curvature. Our direct determination can 
be used to strongly constrain an 0{p'^) covariant baryon 
chiral perturbation theory extrapolation of the nucleon 
mass 1221 . The combined fit yields the preliminary value 



a. 



phys _ 
irN 



(38 ± 12) MeV at the physical point. Without 
the direct information on the slope at one point the statis- 
tical and systematic uncertainties would have been much 
larger, see Fig.[7l It would be difficult to significantly re- 
duce this large error, without nucleon mass data at phys- 
ical and, possibly, smaller than physical quark masses. 

We are also able to exclude values y > 0.14 of the y- 
parameter, with a confidence level of 95 %. This means 
that the strangeness contribution to the scalar coupling 
of the nucleon is much smaller than that due to the light 
quark cr-term. To determine the strangeness cr-term us- 
ing the indirect method, i.e. the Feynman-Hcllmann the- 
orem, requires np — 2 + 1 sea quark fiavours. Even then 
the dependence of the nucleon mass on the strange quark 
mass will be very weak and the tiny slope (and its error) 
will be amplified by the mass ratio ms/mud > 25. There- 
fore, for an accurate prediction of /t^ , and in particular 
for a non-vanishing lower bound on its value, an addi- 
tional direct determination at one or a few mass points 
will be crucial. 

We remark that the results presented here have not 
been extrapolated to the continuum limit. Neither has 
the effect of quenching the strange quark been addressed, 
except within the renormalization of the strangeness cr- 
term. 
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